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It is shown that our multi-high-frequency (40 — 200 GHz) resonant 
cavity technique yields distortion-free high field EPR spectra for sin- 
gle crystal samples of the uniaxial and biaxial spin S = 10 single 
molecule magnets (SMMs) [Mni 2 Oi2(CH3COO)i 6 (H20)4]-2CH 3 COOH-4H 2 
and [Fe802(OH)i2(tacn)6]Brg • 9H2O. The observed lineshapes exhibit a pro- 
nounced dependence on temperature, magnetic field, and the spin quantum 
numbers (M5 values) associated with the levels involved in the transitions. 
Measurements at many frequencies allow us to separate various contributions 
to the EPR linewidths, including significant D— strain, g— strain and broad- 
ening due to the random dipolar fields of neighboring molecules. We also 
identify asymmetry in some of the EPR lineshapes for Fes, an d a previously 
unobserved fine structure to some of the EPR lines for both the Fes an d 
Mni2 systems. These findings prove relevant to the mechanism of quantum 
tunneling of magnetization in these SMMs. 
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PACS numbers: 76.30.-v, 75.50.Xx, 75.45.+j 
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I. Introduction 



Examples of "single molecule magnets" (SMMs) which have attracted consider- 
able recent interest include the uniaxial [Mni 2 Oi 2 (CH3COO)i 6 (H 2 0)4]-2CH 3 COOH4H 2 
(Mn 12 — acetate, or Mn 12 — ac) system, and the biaxial [Fe 8 2 (OH) 12 (tacn) 6 ]Br 8 • 9H 2 
(Fe 8 Br) system |Tjg|||,§§0,§,|rnj . Both of these SMMs exhibit spin S = 10 ground 
states, with a dominant easy-axis type anisotropy, DS% (where D is negative — see section 
II), which provides a potential barrier (height ~ I-DIS" 2 ) against the reversal of individual 
spins from Ms = +10 to —10 (see Fig. 1). Recent studies have shown that both systems 
possess an ability for their magnetic moments to quantum tunnel through their respective 
anisotropy barriers (see Fig. 1), and that this process dominates the magnetization dynamics 
at low temperatures (A; B T< \D\S 2 ) ^^^J^^^jJ^ . 

Because of its significance as a fundamentally novel phenomenon, as well as its potential 
for applications in nano-technology, magnetic quantum tunneling (MQT) in SMMs has been 
extensively studied lately (see e.g. P, p+| , p^2| , p!3| , p^ , p^| . ^6| , p!7| , p!8| ) , yet its underlying mechanism 
is still not fully understood, especially at a quantitative level. One possible reason is the 
lack of precise information on the disposition of the individual (electronic) spin energy levels. 
This is due to the fact that much of the available data on the spin properties of these systems 
have been obtained through DC and AC magnetization studies 0,[l5|,|l6|,|T9| , specific heat 
20fl , neutron scattering pT| , |22|p3| , proton NMR |24]j25|1 , and muon spin relaxation |24j], which 



cannot easily or directly probe the properties of specific spin quantum levels. 

An important breakthrough in the spectroscopy of SMMs has been the development of 
a sensitive cavity perturbation technique for carrying out multi-high-frequency (40 to > 200 



GHz) EPR measurements on single crystals p6| , p7| , p8| , p9| , |30| . In particular, this technique 
leads to at least a three orders of magnitude improvement in detection sensitivity JZB) relative 
to conventional single-pass multi-high- frequency EPR techniques [P, pT] , P^ , P^ , P^ , P5| , P^1 , which 
opens up many new avenues for probing SMMs. To date, our single crystal investigations 
have focused primarily on the energy dependence of the positions in magnetic field of EPR 
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absorptions corresponding to transitions between different spin quantum levels [26| , |27 , |28 , 29 . 
In this article, we turn our attention to the magnetic field, frequency, temperature and spin 
quantum level (Ms) dependence of the EPR linewidths and, in particular, the line shapes 
obtained from measurements on FesBr and Mni2— ac. We note that a brief account has been 
presented elsewhere |J7| . However, this earlier report emphasizes the theoretical foundations 
of the linewidth analysis, and does not consider lineshapes. 

Documented lack of resolution, sensitivity and/or various experimental artifacts have 
been known to affect most high frequency EPR investigations of FesBr and Mni2— ac to 
date. It is commonly stated that observed EPR lineshapes cannot be analyzed because 
they contain uncontrollable distortions due to instrumental artifacts resulting from standing 
waves in the sample probe, which cause a mixing of the dissipative and reactive sample 
response [p2[,|33|,|38f . This is most pronounced in the case of the high magnetic field studies of 



Barra et al. [p| j3^ , |33|j34| ; asymmetric lineshapes are generally observed, making it extremely 
difficult to evaluate the true resonance positions. Indeed, even for aligned polycrystals, 
significant discrepancies can be noted between the observed and calculated EPR transitions, 
with respect to both the peak positions and intensities, as well as the line shapes and 
widths [p3| . As an example, it has been reported in that the linewidths in Fe 8 Br 
approximately scale as Mj, whereas our earlier studies |29|j37|| have shown a dominant linear 
dependence of the linewidths on M$, both for FegBr and Mn 12 — ac. This fact, along with 
the asymmetric lineshapes observed in high field EPR studies, may account for differences 
in the spin Hamiltonian parameters deduced by single-pass (i.e. cavity-less) high-field EPR, 



as compared to other techniques [16,26,28|,|29||. The zero field coherent time-domain THz 



technique employed by Parks et al. | 36[| , is limited in energy resolution (resolution is a fifth 
of the broadest EPR line width), as well as sensitivity. Consequently, it is necessary to study 
a large sample made up from unaligned microcrystals, and it is only possible to observe a 
few transitions. Furthermore, experimental complications only permit a linewidth analysis, 
i.e. no lineshape analysis is possible. To date, the most precise SMM EPR lineshape 
analyses have been reported by Mukhin et al., for both FesBr and Mni2— ac fl35| , using a 



zero-field coherent source technique. Gaussian EPR lines are observed, which are attributed 
to random dipolar fields in the sample. However, these studies are also limited in sensitivity. 
Consequently, only a few transitions are observed, even for very large polycrystalline samples. 
Finally, conventional X-band EPR is limited by the fact that only a few transitions close to 
the top of the barrier are accessible. Furthermore, little or no information concerning the 
spin-spin (T2) and spin-lattice (Ti) relaxation times has been possible in X-band p8 |. 



As we will show in this article, measurements on single crystals are essential, as are 
several other unique aspects of our technique. For example, contrary to recent assertions 
by Blinc et al. ||38|| , we are able to obviate essentially all instrumental artifacts that could 
influence/distort EPR line shapes. In addition, our ability to measure at many frequencies, 
temperatures and orientations enables us to distinguish between several contributions to 
the EPR line widths and shapes. These include: inter-SMM and hyperfine dipolar fields; 
distributions in crystal field parameters; and spin-lattice interactions. Indeed, each of these 



effects have been discussed in MQT literature, e.g. dipolar fields |3J|, distributed crystal 
field parameters caused by dislocations PD| , thermally assisted tunneling ]I3, ffl] , f£2| ) the 



Landau-Zener effect P3||Pi| , fi5| , and hyperfine interactions j|!|[|?J]. Based on our recent 
theoretical work [37], we find that a distribution in the uniaxial crystal field parameter, so 
called D— strain, contributes significantly to the EPR linewidths for the FesBr and Mni2 — ac 
systems. In addition, we observe an appreciable inter-SMM dipolar broadening for FegBr, 
and g— strain in Mn 12 — ac. The current work focuses on the EPR lineshapes. This is 
particularly important in view of recent theoretical work j40j which considers local distortions 
in the SMM site symmetries, and their role in the MQT phenomenon. 

The article is organized as follows: in the next section (II), we discuss the S = 10 
effective spin Hamiltonian in the context of high-frequency EPR and the MQT phenomenon; 
in section III we describe our experimental procedure; in section IV, we separately present 
results for FesBr (IV. A) and Mni2— ac (IV. B); in section V, we discuss the significance of 
our results, and we summarize and conclude our findings in section VI. 
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II. The effective spin Hamiltonian 

The effective spin Hamiltonian for the Mn 12 — ac and FegBr systems has the form: 

H = DS 2 Z + E(S 2 X - SJ) + fMBB-i-S + O4 + H' ; (1) 

where S is the vector spin operator and S z , S x and S y are the spin projection operators 
along the easy, hard and intermediate axes respectively; D is the uniaxial spin-spin coupling 
parameter and the rhombic E (< D/3) term characterizes the magnetic anisotropy in the 
plane perpendicular to the easy axis; g is the Lande g-tensor and B is the applied magnetic 
field vector; finally, O4 denotes weaker higher order terms of order four in the spin operators, 
and H' represents additional perturbations such as those which lead to EPR line broadening, 
as discussed in the introduction. D is negative for both Mn 12 — ac (= —0.457 cm -1 ||22|| ) and 
FegBr (= —0.203 cnr 1 ||23||). The rhombic term E is zero for Mni2— ac, due to its strictly 
axial (S4) site symmetry, while E = 0.032 cm -1 for Fe 8 Br ||23|| . For a cylindrically symmetric 
system, for which DS 2 would be the only non-zero term in Eq. (1) in zero applied field, 
the energy eigenstates may be labeled by the quantum number M s (-S < M s < S), which 
represents the projection of S onto the easy axis; the energy eigenvalues are then given by 
the expression e = DM^. This results in the energy barrier separating doubly degenerate 
(Ms = ±z, i = integer) spin "up" and "down" states (see Fig. 1). 

In order for quantum tunneling to be possible, there must exist finite terms in Eq. (1) 
that break the cylindrical symmetry (i.e. do not commute with S z ) and mix states in either 
potential well, thereby lifting the Ms = ±z degeneracies. Such terms include: the rhombic 
E(S 2 — S 2 ) term; higher order single-ion anisotropy terms in the transverse spin operators 
(e.g. S x and Sy) [|],|16|,^2j ; a transverse externally applied magnetic field; transverse internal 
fields due to neighboring SMMs or nuclei (contained in H' p9|^6| ); or local distortions in 
the crystal field symmetry caused imperfections in the crystal (also contained in H' pi) 



— these may include dislocations, or a partial loss or disorder among the ligand molecules. 
Pure quantum tunneling is observed between the unperturbed Ms = ±10 degenerate ground 
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states in FesBr below about 350 mK (see Fig. 1), and the observed relaxation has been 



explained in terms of the E(S% — Sy) term in Eq. (1) p£B |. However, one has to additionally 
consider the respective roles of inter-SMM dipolar fields and nuclear hyperfine fields [H 1 
in Eq. (1)] in order to fully account for the observed tunneling in FesBr [Hj|,£6|j. For the 
Mn 12 — ac system, the situation is even more complex: pure ground state quantum tunneling 
has not been observed in zero applied magnetic field, and there is little consensus as to the 
origin of the apparent (weak) transverse anisotropy which causes thermally assisted quantum 
tunneling among levels near the top of the barrier (see Fig. 1). The leading theories include: 
fourth order single ion transverse anisotropy terms such as S% + Sy, which are allowed under 
the axial (S4) site symmetry [|3| |T6"l , p2f ; and local distortions in the axial crystal fields caused, 



e.g. by dislocations [^0 



A tunneling matrix element (term in Eq. 1 that does not commute with S z ) gives rise 
to eigenstates that are constructed from symmetric and antisymmetric combinations of un- 
perturbed (Ms = ±2) states on either side of the barrier, and to a lifting of the degeneracy 
(tunnel splitting) between these symmetric and antisymmetric states. Thus, the tunneling 
matrix elements modify the energy eigenvalues and should, therefore, be detectable in an 
EPR experiment. However, for the systems of interest, these effects are small. Consequently, 
great precision is required. Although inelastic neutron scattering (INS) measurements are 
capable of providing similar information, EPR offers many advantages in terms of resolution, 
sensitivity (amount of sample required) and, in particular, the fact that one can systemati- 
cally control the level mixing by means of an externally applied magnetic field. In addition, 
the higher sensitivity of the EPR technique assures that all of this may be achieved for a 
single crystal. Several EPR (and INS) experiments have indicated the presence of significant 



fourth order single- ion transverse anisotropies for both Mn 12 — ac [ 16, 26, 32, 35 1 and Fe§Br 



33| , |35[| . Nevertheless, fourth order anisotropy cannot fully explain the observed spacing in 



magnetic field of magnetization steps observed in low temperature hysteresis experiments 
, and there remains considerable disagreement as to the magnitude of the fourth order 



terms [p^,[^,^,[^,0.[||,|5T|]. Consequently, higher resolution EPR measurements are neces- 



sary in order to refine estimates of the Hamiltonian parameters, and to determine whether 
EPR lineshapes contain additional information concerning the MQT phenomenon. 



III. Experimental 

The (2S + 1)— fold quantum energy level structure associated with a large molecular spin S 
necessitates spectroscopies spanning a wide frequency range. Furthermore, large zero-field 
level splittings, due to the significant crystalline anisotropy (large \D\) and large total spin S, 
demand the use of frequencies and magnetic fields considerably higher (50 GHz to 1 THz, and 
up to 10 tesla respectively) than those typically used by the majority of EPR spectroscopists. 
The high degree of sensitivity required for single crystal measurements is achieved using a 
resonant cavity perturbation technique in combination with a broad-band Millimeter-wave 
Vector Network Analyzer (MVNA) exhibiting an exceptionally good signal-to-noise ratio; 
a detailed description of this instrumentation can be found in ref. p8|j5^| . The MVNA is 
a phase sensitive, fully sweepable (8 to 350 GHz |53[), superheterodyne source/detection 
system. Several sample probes couple the network analyzer to a range of high sensitivity 
cavities (Q— factors of up to 25,000) situated within the bore of a 7 tesla superconducting 
magnet. The MVNA/cavity combination has been shown to exhibit a sensitivity of at least 
10 9 spins-G _1 -s _1 , which is comparable with the best narrow-band EPR spectrometers. 
This, coupled with newly acquired sources and a split-pair magnet, allow single crystal 
measurements at any frequency in the range from 8 to 250 GHz, at temperatures down to 
1.5 K (±0.01K), and for any geometrical combination of DC and AC field orientations up 
to 7 T (up to 45 T at the National High Magnetic Field Laboratory). 

The use of a narrow band cavity offers many important advantages over non-resonant 
methods (see e.g. refs B, |5I] , |5^ , |53| , |53}| ). Careful consideration concerning the coupling of 
radiation to and from the cavity (via waveguide), combined with the ability to study very 
small samples, eliminates problems associated with standing waves in the sample probe 
|52|]. This, in turn, eliminates a mixing of the dissipative and reactive responses of the 



sample under investigation and, when combined with a vector detection scheme, enables 
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faithful extraction of the true EPR lineshapes (both the real and imaginary components), 
free from instrumental artifacts. One other notable feature of the superheterodyne scheme 
is its detection rate of 34 MHz; thus, there is no need for field modulation, unlike other 
high field EPR methods P J31|J3^J33[| . Consequently, any raw data displayed in this paper 
constitutes pure absorption. Finally, the use of a cavity enables positioning of a single 
crystal sample into a well defined electromagnetic field environment, i.e. the orientations 
of the DC and AC magnetic fields relative to the sample's crystallographic axes is precisely 
known, and the DC field is very homogeneous over the tiny volume (<C 1mm 3 ) of the sample. 
In this way, additional contributions to the EPR lineshapes are avoided due to, e.g. field 
inhomgeneities over the sample volume, or slight mis-alignments of crystallites, as may be 
the case for measurements on aligned powders. 

Single Mn 12 — ac and FegBr crystals were grown using literature methods 0,0. All mea- 
surements were performed in the standard EPR configuration with the AC excitation field 
transverse to the DC field. The FegBr crystal used for the majority of the studies reported 
in the following section was a thin rhombic (acute angle ~ 60°) platelet with edges of length 
~ 0.7 mm, and thickness ~ 0.17 mm. The Mn 12 — ac crystal was needle shaped, with approx- 
imate dimensions 1 x 0.15 x 0.05 mm 3 . Orientation of Mn 12 — ac crystal was straightforward, 
due to the fact that its easy axis is defined by the needle axis of the sample. Alignment of 
the FegBr crystal necessitated measurements at many angles |[29|| . All of the data presented 
in this paper were obtained with the DC magnetic field aligned approximately parallel to the 
samples' easy axes (hard axis data are presented elsewhere p6| . |29| ). Finally, temperature 
control, in the range 2 to 50 K, was achieved using a variable flow cryostat situated within 
the bore of the superconducting magnet. 

IV. Results 

A. FegBr 

Figure 2 shows raw temperature dependent data obtained for the FegBr SMM at a frequency 
of 116.931 GHz; the field is applied along the easy axis to within an accuracy of about 1°. 
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A series of more-or-less evenly spaced resonant absorptions are observed which become 
successively weaker with increasing field. A detailed analysis of the frequency dependence 
of the EPR line positions is published elsewhere ||29|| . Here we focus instead on the line 
widths and shapes. In this frequency and field range, and for this orientation, the levels 
involved in the transitions are well separated and do not cross. Thus, M5 is a good quantum 
number, and we have labeled the resonances in Fig. 2 according to the levels involved in 
the transitions. 

It is apparent from Fig. 2 that the observed transitions exhibit dramatic variations in 
linewidth, both as a function of Ms, field and temperature. In order to characterize this 
width dependence, it is first necessary to fit the line shapes. Fig. 3 shows several attempts 
to fit various transitions with both Lorentzian and Gaussian functions; all of the data were 
obtained from a single field sweep at a frequency of 89.035 GHz and a temperature of 10 K, 
with the field close to alignment with the easy axis (8° ± 1°). At the lowest fields, where the 
resonances are broad (comparable to spacings) and therefore overlap, fits involving multiple 
Gaussian or Lorentzian functions are necessary (Fig. 3). As can be seen from this figure, 
the Gaussian fit is noticeably better than the Lorentzian one, suggesting that some form of 
inhomogeneous broadening mechanism is at play. Although this has been noted by other 
authors these earlier measurements were performed on polycrystalline powders, which 
were pressed into large pellets. Consequently, the microcrystals in the pellet will have been 
subjected to varying stresses that could be responsible for the broad Gaussian lines. Here, 
we demonstrate that this is an intrinsic feature of the Feg SMMs (this is also the case for 
Mn 12 — see following section). 

Fig. 4 shows Lorentzian and Gaussian fits to transitions excited from successively higher 
Ms levels (smaller absolute value of M5) from the same field sweep as the data in Fig. 3. 
The first thing to note is the dramatic line narrowing with decreasing absolute value of M^. 
Also apparent (inspite of the reduced signal-to-noise) is a slight change in lineshape, from 
approximately Gaussian to a more Lorentzian form, as the lines become narrower. The 
significance of this finding will be discussed further below. 
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Figure 5 shows a compilation of Gaussian widths, versus Ms (the level from which the 
transition was excited), measured at many different frequencies; all of the data were obtained 
at a temperature of 10 K, and with the field applied at an angle 9 = 8° ± 1° away from 
the easy axis. The pronounced (almost linear) increase in EPR linewidth with Ms is due 
to D— strain [37]. Since the dominant contribution to the spin Hamiltonian (Eq. 1) is 



proportional to M|, energy differences {i.e. EPR transition frequencies) scale as Ms; hence 
D— strain produces a linear width (Aenergy) dependence on Ms, which projects onto a 
width in field (AB) that also scales linearly with Ms- The rounding close to Ms = is due 
to a convolution of the intrinsic lifetime (spin-lattice) broadening and the Ms dependent 
(D— strain) contribution. The former is expected to produce a Lorentzian lineshape, thus 
accounting for the apparent crossover from Gaussian to Lorentzian EPR lines, as observed 
for the sharper transitions in Fig. 4. The slight narrowing of the EPR widths in Fig. 5 as a 
function of increasing frequency, and the weak asymmetry about Ms = —1/2, is due to the 
fact that higher frequency (higher Ms) transitions are observed at higher magnetic fields, 
where the inter-SMM dipolar broadening is weaker (see discussion below). For Fe§Br, no 
measurable g— strain is observed. A theoretical analysis of these trends may be found in ref. 



The temperature dependent contribution to the linewidths for FesBr (Fig. 6) is domi- 
nated by random/fluctuating dipolar fields at each SMM site, due to the magnetic moments 
of the neighboring SMMs. At 116 GHz, the Ms = —10 to —9 transition occurs very close 
to zero field (~ 0.1 T). Under zero field cooling, the spin- up (Ms = —10) and spin-down 
(Ms = +10) states will be equally populated. Hence, the random distribution of spin-up 
and spin-down molecules results in a maximal contribution to the Gaussian dipolar EPR line 
broadening, as evidenced by the extremely large temperature dependent contribution to the 
Ms = —10 to —9 transition in Fig. 6 as T— ► 0. Upon raising the temperature, higher levels 
in the scheme (|Ms| < 10) become populated and the individual moments begin to fluctuate 
more rapidly, thus leading to a suppression of the dipolar broadening. The same qualita- 
tive picture may be applied to the other transitions (Ms = —9 to —8 etc.), provided that 
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ksT exceeds the Zeeman splitting of the ground state. For fc#T less than the ground state 
Zeeman splitting, the system polarizes and the dipole fluctuations are suppressed, leading 
to an EPR line narrowing, as confirmed for all but the Ms = — 10 to —9 transition in Fig. 
6. The crossover occurs when the Ms = +10 state starts to depopulate, i.e. approximately 
when fcgT = gSfig3. The inset to Fig. 6 shows the crossover temperature (T C r), where the 
slope of the data in the main part of the figure changes sign, versus the average magnetic 
field strength at which each particular transition is observed (see Fig. 2). It is found that 
Ten indeed scales linearly with B. Furthermore, the slope of the fit to the data in the inset 
to Fig. 6 is 14.4 ± 1 K-T -1 = gS^s/kg, corresponding very well to S — 10 and g ~ 2 to 
within the experimental error, thus supporting the picture of dipolar broadening. T C r for 
the 116 GHz Ms = — 10 to —9 transition in Fig. 6 is below the minimum temperature of our 
current spectrometer (T min ~ 1.2 K). Consequently, we were unable to verify whether this 
line eventually narrows when the spin system polarizes. However, higher frequency /field 
data, for which Ten should exceed the minimum temperature of our spectrometer, exhibit a 
similar line broadening for the Ms = — 10 to —9 transition, at all temperatures investigated 
(see Fig. 8 below). At present, the origin of this effect is not fully understood. 

In Fig. 7, we show the temperature dependent shifts of the EPR line positions obtained 
from Fig. 2. When compared to Fig. 6, it is noticeable that there is a similar crossover 
in the temperature dependence for transitions from levels with \Mg\ < 10. Therefore, this 
low-temperature limiting behavior may also be attributed to the polarization of the spin 
system. The sign of the shift indicates that the field at each SMM site, due to neighboring 
SMMs, opposes the applied field, i.e. it is opposite to the polarization direction. This is 
to be expected for the FesBr structure, in which the Fes molecules reside on the corners 



of a slightly distorted cube From the magnitudes of the |M,g| < 10 line shifts, we 

estimate dipolar fields at each SMM site of a few hundred Gauss (~ 300 G). The origin of 
the continued shift of the Ms = —10 to —9 transition to lower magnetic fields is, at present, 
not well understood. 

A comprehensive theoretical analysis of this linewidth data has been published elsewhere 
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37J . From such an analysis for the sample used in obtaining the data in Fig. 2, we esti- 
mate a distribution in D of <td = 0.01D, negligible (/—strain, and a dipolar contribution 
to the broadening which assumes a nearest neighbor SMM separation of 12 A |[55|| , which 
is in agreement with the accepted range of near neighbor distances of between 10.5 — 14A 
Within acceptable ranges, no other combination of parameter values is able to produce the 
same quality of fit simultaneously to all of the frequencies and temperatures studied for this 
particular sample, even though other parameter values may produce good fits for a partic- 
ular measurement. This emphasizes the power of our multi- frequency single crystal EPR 
technique. 

Very recent studies on a different sample from the one used in the measurements described 



above (similar shape and size, and obtained from the same synthesis [[56]]), reveal unusual 
lineshape characteristics, particularly at low temperatures and for transitions involving levels 
close to the bottom of the potential wells (e.g. Ms = —10 — > —9 and = — 9 — > —8,..). 
This is illustrated in Fig. 8, for data obtained at a frequency of 145.937 GHz, with the field 
aligned precisely along the easy axis. At the somewhat elevated temperature of 10 K (Fig. 
8a), the Ms = — 10 — > —9 transition shows appreciable asymmetry: this resonance has a 
Lorentzian (dotted curve) tail on the low field side, and closer to a Gaussian (dashed curve) 
tail on the high field side; furthermore, the data lie somewhat to the right of the fits on the 
high field side. Subsequent resonances (M s = — 9 — ► —8, M s = — 8 — > —7, etc.) appear to be 
less asymmetric, and the Gaussian fit improves. Recent theoretical studies have considered, 
in some detail, the effects of different types of defect structure on EPR lineshapes, and 
how these defects may influence the MQT phenomenon in SMMs (see section V) |40| , |57| . 
Both Lorentzian and Gaussian broadening are predicted [f40"| , as well as asymmetry |57j (see 
discussion in section V). 

Examination of the = —10 — > —9 transition at lower temperatures (Fig. 8b) reveals 
a fine structure to the line, which is strongly temperature dependent, i.e. the curves are 
not smooth, but exhibit a series of distinct distortions. We suspect that this fine structure 
and the asymmetry is connected with crystal defects, since it is not observed for all samples 
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56fl . A direct comparison between the data in Fig. 8 and data for the first sample (Figs 
3 and 4) is complicated by the fact that the linewidths are substantially broader for this 
second sample (increased D— strain). We note that many SMM compounds are known to 
lose solvent molecules pEIJjBfl , which would lead to a degredation of sample quality over 
time and may, in fact, provide a natural explanation for the D— strain. Nevertheless, the 
Ms = —10 — > —9 transition in Fig. 2 (lowest field resonance) also exhibits fine structure 
at the lowest temperatures, and fine structure is also apparent in data published by Barra 
et al. []33 |. It is possible that the fine structure is indicative of a finite number of SMM 



site symmetries, as opposed to a smooth distribution. The dislocation theory of Garanin 
and Chudnovsky |4(J leads to long range strains in the crystal, thus giving rise to a smooth 



distribution in site symmetries. Therefore, the fine structure we observe may signify a finite 
number of configurations involving Fe§ molecules coordinated by different numbers or ar- 
rangements of ligand molecules. It should be stressed that the data presented in Fig. 8 truly 
represent the intrinsic EPR lineshapes, i.e. any instrumental effects would be systematic 
and, therefore, influence all of the lineshapes equivalently. Furthermore, the noise level for 
these measurements is well below the resolution of the figure. Clearly, additional studies on 
a range of samples prepared in different ways, i. e. different sample qualities, will be required 
to fully resolve these issues — such studies are currently in progress. 

B. Mni 2 -ac 

Figure 9 shows raw data obtained for Mn 12 — ac at a frequency of 189.123 GHz, and at 
temperatures between 10K and 35K; the field is applied parallel to the easy axis to within 
0.2°. The first point to note is the diminished signal-to-noise for this data, when compared 
to Fig. 2. There are several reasons for this: first, the dynamic range of our spectrometer 
is reduced at this higher frequency of 189 GHz; second, the volume of the sample used in 
this investigation was somewhat smaller than that of the FesBr sample used for the studies 
in the previous section; finally, due to the increased energy scale of the Mn 12 — ac system 
(larger \D\), the observed transitions correspond to levels considerably higher in energy 
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above the ground state than would be the case for FegBr at comparable frequencies and 
fields. Therefore, this constitutes a much more challenging experiment. Nevertheless, clear 
resonances are observed, which also have Gaussian lineshapes (solid curves), and once again 
exhibit a pronounced linewidth dependence on Ms (the level from which the transition was 
excited). For comparison, the inset to Fig. 9 shows attempts to fit the first two transitions 
with Gaussian and Lorentzian functions; the Gaussian fit is noticeably better. 

Figure 10 plots the Ms dependence of the Gaussian EPR linewidths for measurements 
made at several different frequencies. The linear Ms contribution to the widths once again 
indicates a significant distribution in D (o£> ~ 0.02D). However, the dipolar contribution to 
the linewidths is considerably weaker for Mn 12 — ac than for FegBr. This is initially appar- 
ent from the weak frequency dependence observed amongst the widths for transitions with 
| Ms | < 5, and the fact that the frequency dependence observed for higher |M^| transitions 
is opposite to what was observed for Fes in Fig. 5, i.e. the higher frequency transitions 
are broader. The origin of this effect has to do with g-strain, and will be discussed further 
below. 

In Fig. 11, we plot the temperature dependence of the linewidths observed at a fre- 
quency of 189.123 GHz. A similar trend can be observed for the low |M^| transitions for 
both Mni2— ac and FegBr (decreasing width with decreasing temperature). However, the 
pronounced line narrowing observed for the higher |Ms| transitions at lower temperatures, 
in FegBr, is absent in Mn 12 — ac. This provides further evidence that the dipolar contribution 
to the linewidths is considerably weaker in Mn 12 — ac. In contrast, the Ms dependence of the 
widths is stronger for Mn 12 — ac, due to the increased D— strain (a D ps 0.02D, as opposed 
to od ~ 0.01/} for FegBr). The weaker dipolar broadening is consistent with an increased 
inter-SMM separation for Mn 12 — ac [0. 

The pronounced frequency dependence of the larger |Ms| (> 5) transitions, as seen in 
Fig. 10, is attributable to g— strain. Each transition shifts to higher field with increas- 
ing frequency, thereby amplifying the effect of a distribution in g — hence the increased 
linewidths with increasing frequency. Additional support for g— strain comes from the mea- 

15 



sured anisotropy in the g-tensor, both from our previous work J26] and that of others |3T| , |32 • 



Using accepted values for the inter-SMM separations, it is possible to simulate all data ob- 
tained for a particular sample with a single pair of (sample dependent) D— and g— strain 
parameters; od = 0.02D and a g = 0.008g are typical for many of the Mni2— ac samples we 
have studied P7| . 

Finally, in Fig. 12, we show magnifications of several of the resonances from Fig. 9. 
Once again, it is apparent that the EPR lines exhibit a fine structure. In particular, the 
Ms = —6 to —5 and Ms = —5 to —4 transitions show a clear double peak structure; the best 
attempts at single Gaussian fits (dotted curves) are included in the figure for comparison. 
Similar behavior was observed at other frequencies. Once again, these findings hint at 
the possibility that there may exist a finite number of configurations involving Mn 12 — ac 
molecules coordinated by different numbers or arrangements of ligand molecules. We note 
that the presence of various isomeric forms of Mni2 , whereby the arrangements of the ligands 
differ slightly, is well established both in the acetate |59[ and in Mn 12 SMMs involving other 
ligands 



V. Discussion 

The unambiguous demonstration of D— strain in FesBr and Mni2— ac indicates the presence 
of significant sample inhomogeneities or imperfections, i.e. the molecules are subjected to 
varying anisotropy. A natural explanation for this could be a loss of solvent molecules, or 
disorder among the ligands; indeed a slight degredation of sample quality over time has been 



noted, and it is known that the crystals do lose solvent |58| . Another possibility has recently 
been proposed by Garanin and Chudnovsky [^0[ , involving dislocations in the crystals which 
act to locally distort the site symmetries of individual SMMs. What is more, Garanin and 
Chudnovsky have claimed that the QTM phenomenon in SMMs finds a natural explanation 
in terms of these dislocations. As mentioned in the introduction, QTM requires symmetry 
breaking terms in the Hamiltonian (Eq. 1). Local rotations of the magnetic anisotropy axes 
due to dislocations have two effects: they generate transverse second order crystal field terms, 
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of the form E(S^, — Sy), which facilitate QTM between Ms states differing by an even integer; 
and they result in an effective local transverse magnetic field (under application of field || z) 
that unfreezes odd tunneling resonances, thereby solving the conundrum as to why both 
even and odd tunneling resonances are always observed in FegBr and Mn 12 — ac. It has also 
been shown that long-range strains produced by dislocations result in broad distributions of 
relaxation times, something which is also observed experimentally in Mni2— ac [(5J. Indeed, 
assuming realistic concentrations of dislocations, Garanin and Chudnovsky [[KJ have been 
able to account for the experimentally observed tunneling in Mn 12 — ac [fUSfl , and predict a 
distribution in the uniaxial spin-spin coupling parameter od w 0.027.D, which is of the same 
order as we have found from linewidth analyses (<7d ~ 0.01 — 0.02D). It should be pointed 
out, however, that any disorder amongst the ligand molecules, as well as any loss of solvent 
molecules, may be expected to produce many of the same effects as those proposed by 
Garanin and Chudnovsky f40fl , i.e. disorder amongst the peripheral molecules surrounding 
the SMMs will result in local variations in site symmetry. 

Careful theoretical consideration of different forms of dislocations (edge, screw, corre- 
lated, uncorrelated, etc..) indicate that the resultant D— strain can produce both Gaussian 
and Lorentzian EPR lineshapes ||40|| . Furthermore, recent work by Park et al. |57||, predicts 
asymmetry in some of the EPR lineshapes; one may think of this as arising due to the fact 
that a dislocation causes an asymmetric distribution in D, i.e. it can only reduce D, not 
enhance it. Clearly, certain features of the results presented here, particularly in Fig. 8a, 
may confirm some of these predictions. So far, our attempts to simulate data have been 
limited to D— strain, g— strain and dipolar broadening. The D— and g— strains should be 
correlated, although we have yet to verify this. Simulated E— strain, as well as strain in 
the fourth order anisotropy, did not affect our simulations to within sensible ranges of these 
parameters. We have also ignored hyperfine effects, since these are irrelevant in the case 
of the FegBr ( 56 Fe has no nuclear spin), and exchange narrowing in Mn 12 — ac is believed to 
diminish the hyperfine splitting |37fl. At present, most of our data have been obtained for 
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just a few samples of each of FesBr and M1112— ac. Further work is in progress on a range of 
different SMMs, with the goal of resolving some of these issues, e.g. whether solvent loss or 
ligand disorder is the main reason for the observed D— strain, or whether EPR lineshapes 
correlate with various theoretical works in progress dealing with dislocation mediated QTM 



Clearly, if the QTM mechanism is mediated by disorder, then single crystal EPR offers 
the most direct means of providing quantitative information with which to test relevant the- 
oretical predictions. However, it is essential that EPR lineshapes obtained from experiments 
are free from the experimental artifacts which have been known to hamper high frequency 
EPR measurements. Indeed, given the strong Mg dependence of the EPR line shapes and 
widths, one has to question recent attempts to determine the crystal field parameters to 
fourth order from high frequency EPR. Inspection of the fits to data for FegBr and Mn 12 — ac 
by Barra et al. P, p2| , p3| , p4[ , reveal considerable discrepancies in the positions of low |Ms| 
transitions. Because of the smaller linewidths associated with the low |Ms| transitions, one 
can be far more confident as to the experimentally determined resonance field from such 
transitions, than one can be for the considerably broader high |Mg| resonances. Field mod- 
ulation techniques compound this problem, as do studies on aligned polycrystals which give 
rise to asymmetric distributions in the orientations of the individual crystallites. Indeed, 
our most recent investigations indicate that the high frequency EPR technique described in 
this article (section III and [ 28|,52| ) yields crystal field parameters that are in much closer 



agreement with inelastic neutron scattering measurements than conventional high field EPR 



studies 61 



VI. Summary and conclusions 

We have used a multi-high frequency cavity perturbation technique to measure EPR 
lineshapes for oriented single crystal samples of the FesBr and Mni2— ac SMMs. The use of 
a narrow band cavity offers many important advantages over non-resonant methods: careful 
consideration concerning the coupling of radiation to and from the cavity (via waveguide), 
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combined with the ability to study very small samples, eliminates mixing of the dissipative 
and reactive responses of the sample under investigation, and enables faithful extraction of 
the true EPR lineshapes. Our ability to measure at many frequencies, temperatures and 
orientations further allows us to distinguish between several contributions to the EPR line 
widths and shapes. These include: inter-SMM and hyperfine dipolar fields; distributions 
in crystal field parameters; and spin-lattice interactions. Each of these effects have been 
discussed in the literature in terms playing a possible role in the MQT phenomenon. In 
addition to significant D— strain, which is indicative of a distribution in the SMM site sym- 
metries, we have detected a previously unobserved fine structure to some of the resonances 
for both Fe 8 Br and Mni 2 -ac SMMs. 
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Figure captions 



magnetic field, with a dominant easy axis crystalline anisotropy DS^, with D < 0; the M 



Fig. 1 a) Schematic of the double potential well for a spin S — 10 system in zero applied 

s = 

±10 "spin-up" and "spin-down" states lie lowest in energy. Pure QTM, directly between 
the unperturbed Ms = ±10 states (as indicated), is possible if there exist terms in the 
Hamiltonian (Eq. 1) which mix these two levels. In the absence of pure QTM, a thermally 
assisted form of QTM may occur if states higher up in the scheme are mixed, b) Application 
of a magnetic field biases the wells, favoring population of the — Ms states. Whenever levels 
in the two wells come into resonance, direct QTM is possible, and accelerated relaxation 
occurs between the spin-up and the spin-down states, i.e. the relaxation may be considerably 
faster than the classical thermally activated relaxation over the top of the barrier. Resonant 



QTM has been observed in the form of sharp steps in the hysteresis loops of SMMs |jl6 



Fig. 2. Temperature (indicated in the figure) dependent microwave absorption amplitude 
for Fe§Br, as a function of magnetic field, for measurements at 116.931 GHz; the field is 
applied parallel to the easy axis to within an accuracy of 1°. Sharp resonances are observed, 
which correspond to EPR transitions between the 25* + 1 spin states of the system — these 
are labelled in the figure. 

Fig. 3. A close up of several FegBr resonances obtained at 10 K and a frequency of 89.035 
GHz, with the field applied approximately parallel to the easy axis; the levels involved in 
the transitions are indicated in the figure. The solid and dashed curves represent respective 
Gaussian and Lorentzian fits to the data. 

Fig. 4. Further Gaussian and Lorentzian fits to FesBr data (same experimental condi- 
tions as Fig. 3), illustrating the Mg (level from which the transition was excited — indicated 
in the figure ) dependence of the widths. The horizontal scale has been offset and normalized 
to the widths of each resonance — the actual widths are indicated in the figure. 

Fig. 5. A compilation of the Ms (level from which the transition was excited) dependence 

24 



of the Gaussian EPR linewidths obtained for Fe 8 Br at 10 K and many different frequencies 
(indicated in the figure); see text for a discussion of the data. 

Fig. 6. Temperature dependence of the Gaussian EPR linewidths obtained for different 
spin transitions (indicated in the figure); the data were obtained at a frequency of 116.931 
GHz, with the field applied parallel to the easy axis. The inset plots the field dependence 
of the crossover temperature (Tcr) ~ the temperature at which the linewidths reach a 
maximum — for the Ms = —9 to —8, —8 to —7, —7 to —6 and —6 to —5 transitions. 

Fig. 7. Temperature dependence of the EPR line positions obtained for different spin 
transitions (indicated in the figure); the data were obtained at a frequency of 116.931 GHz, 
with the field applied parallel to the easy axis. The line shifts [/x H(T)] are plotted relative 
to the respective line positions at 50 K [/x o H(50 K)]. 

Fig. 8a) A close up of several Fe 8 Br resonances obtained at 10 K and a frequency of 
145.937 GHz, with the field applied exactly parallel to the easy axis; the levels involved in 
the transitions are indicated in the figure. The solid and dashed curves represent respective 
Gaussian and Lorentzian fits to the data, b) Temperature dependence of the M s = —10 to 
—9 transition in (a) illustrating the fine structure to this line. 

Fig. 9. Temperature (indicated in the figure) dependent microwave absorption amplitude 
for Mn 12 — ac, as a function of magnetic field, for measurements at 189.123 GHz; the field is 
applied parallel to the easy axis. Sharp resonances are observed, which correspond to EPR 
transitions between the 2S + 1 spin states (indicated in the figure) of the system. The solid 
curves are Gaussian fits to the data; the inset shows attempts to fit both Gaussian (solid 
curve) and Lorentzian (dashed curve) functions to some of the data in the main part of the 
figure. 

Fig. 10. A compilation of the Ms (level from which the transition was excited) depen- 
dence of the Gaussian EPR linewidths obtained for Mn 12 — ac at 20 K and several different 
frequencies (indicated in the figure); see text for a discussion of the data. 
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Fig. 11. Temperature dependence of the Gaussian EPR linewidths obtained for several 
spin transitions (indicated in the figure) in Mn 12 — ac; the data were obtained at a frequency 
of 189.123 GHz, with the field applied parallel to the easy axis. For comparison, see data 
for Fe 8 Br in Fig. 6. 

Fig. 12. Magnifications of several of the 189.123 GHz traces shown in Fig. 9 (transitions 
and temperatures are indicated in the figure), illustrating a fine structure. 
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FIG. 4. Hill et al. 
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FIG. 8. Hill et al. 
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